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ABSTRACT 

We investigate QCD coherence effects in deep in- 
elastic scattering (DIS) off nuclei and in Drell-Yan 
(DY) dilepton production in proton-nucleus collisions 
within the light-cone color-dipole approach. The 
physical mechanisms underlying the nuclear effects be- 
come very transparent in this approach and are ex- 
plained in some detail. We present numerical calcula- 
tions of nuclear shadowing in DIS and DY and com- 
pare to data. Nuclear effects in the DY transverse 
momentum distribution are calculated as well. The 
dipole approach is the only known way to calculate 
the Cronin effect without introducing additional pa- 
rameters for nuclear targets. 

Keywords: Deep Inelastic Scattering, Drell-Yan Pro- 
cess, QCD, Nuclear Effects, Coherence Effects, Multi- 
ple Scattering. 

1. INTRODUCTION 

The use of nuclei instead of protons in high energy 
scattering experiments, like deep inelastic scattering, 
provides unique possibilities to study the space-time 
development of strongly interacting systems. In ex- 
periments with proton targets the products of the 
scattering process can only be observed in a detec- 
tor which is separated from the reaction point by a 
macroscopic distance. In contrast to this, the nuclear 
medium can serve as a detector located directly at 
the place where the microscopic interaction happens. 
As a consequence, with nuclei one can study coher- 
ence effects in QCD which are not accessible in DIS 
off protons nor in proton-proton scattering. 




FIG. 1: At low XBj and in the target rest frame, the virtual 
photon 7* converts into a qq-palv long before the target. 

At high energies, nuclear scattering is governed by 
coherence effects which are most easily understood in 



the target rest frame. In the rest frame, DIS looks like 
pair creation from a virtual photon, see Fig. 1. Long 
before the target, the virtual photon splits into a qq- 
pair. The lifetime l c of the fluctuation, which is also 
called coherence length, can be estimated with help of 
the uncertainty relation to be of order ~ 1/m^XBj, 
where xsj is Bjorken-x and tun ~ 1 GeV is the nu- 
cleon mass. The coherence length can become much 
greater than the nuclear radius at low XBj- Multi- 
ple scattering within the lifetime of the qq fluctuation 
leads to the pronounced coherence effects observed in 
experiment. 

The most prominent example for a coherent inter- 
action of more than one nucleon is the phenomenon 
of nuclear shadowing, i.e. the suppression of the 
nuclear structure function F£ with respect to the 
proton structure function F% at low XBj < 0.1, 
F^{x Bj , Q 2 )/{AF^{x B j,Q 2 )). Here Q 2 is the virtual- 
ly of the photon. Shadowing in low XBj DIS and at 
high photon virtualities is experimentally well studied 
by NMC [1]. The analogous effect for DY at low x 2 
was discovered by E772 [2]. 

What is the mechanism behind this suppression? 
If the coherence length is very long, as indicated in 
Fig. 1, the gg-dipole undergoes multiple scatterings 
inside the nucleus. The physics of shadowing in DIS 
is most easily understood in a representation, in which 
the pair has a definite transverse size p. As a result of 
color transparency [3, 4, 5], small pairs interact with 
a small cross section, while large pairs interact with 
a large cross section. The term "shadowing" can be 
taken literally in the target rest frame. Large pairs are 
absorbed by the nucleons at the surface which cast a 
shadow on the inner nucleons. The small pairs are not 
shadowed. They have equal chances to interact with 
any of the nucleons in the nucleus. From these simple 
arguments, one can already understand the two nec- 
essary conditions for shadowing. First, the hadronic 
fluctuation of the virtual photon has to interact with 
a large cross section and second, the coherence length 
has to be long enough to allow for multiple scattering. 

While the target rest frame picture of DIS is 
very popular, the light-cone approach of Kopeliovich 
[6, 7, 8], which describes the DY process [9] in the 
target rest frame, is less known. In the light-cone ap- 
proach, DY dilepton production in the rest frame of 
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FIG. 2: A quark (or an antiquark) inside the projectile 
hadron scatters off the target color field and radiates a 
massive photon. The subsequent decay of the 7* into the 
lepton pair is not shown. 



the target appears as bremsstrahlung, see Fig. 2. A 
quark from the projectile scatters off the target and 
radiates a virtual photon. This photon decays into a 
lepton pair. Even though this mechanism looks quite 
different from the usual annihilation picture, it yields 
numerical values very similar to the next-to-leading 
order parton model [11]. Note that while cross sections 
are Lorentz invariant, the partonic interpretation of a 
hard scattering process is not. That is why a process 
that looks like quark-antiquark annihilation in the in- 
finite momentum frame can appear as bremsstrahlung 
in the target rest frame. 

Remarkably, the DY cross section can be expressed 
in terms of the same dipole cross section that appears 
in DIS. In DY, the dipole cross section arises from the 
interference of the two graphs in Fig. 2. On a nuclear 
target, the quark will of course scatter several times. 
The effect of multiple scattering on bremsstrahlung 
is well known in QED as the Landau-Pomcranchuk- 
Migdal (LPM) effect [10]. The LPM effect leads to a 
reduction of the cross section due to destructive inter- 
ferences between the amplitudes for radiation off dif- 
ferent points along the quark's trajectory. In the case 
of DY, this suppression is observed in experiment as 
nuclear shadowing [2]. 



2. SHADOWING IN DIS 

First, we discuss the case of shadowing in deep in- 
elastic scattering off nuclei. Like shadowing in hadron- 
nuclcus collisions [12], also shadowing in DIS is inti- 
mately related to diffraction. The close connection be- 
tween shadowing and diffraction becomes most trans- 
parent in the formula derived by Karmanov and Kon- 
dratyuk [13]. In the double scattering approximation, 
the shadowing correction can be related to the diffrac- 
tion dissociation spectrum, integrated over the mass, 
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is the formfactor of the nucleus, which depends on the 
coherence length 



lc = 
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v is the energy of the 7* in the target rest frame and 
Mx is the mass of the diffractively excited state. The 
coherence length can be estimated from the uncer- 
tainty relation and is the lifetime of the diffractively 
excited state. If l c —> 0, the shadowing correction in 
Eq. (1) vanishes and one is left with the single scat- 
tering term Aa 1 p . 

Note that Eq. (1) is valid only in double scatter- 
ing approximation. For heavy nuclei, however, higher 
order scattering terms will become important. These 
can be calculated, if one knows the eigenstates of the 
interaction. Fortunately, the eigenstates of the T ma- 
trix (restricted to diffractive processes) were identified 
a long time ago in QCD [3, 14] as partonic configura- 
tions with fixed transverse separations in impact pa- 
rameter space. For DIS, the lowest eigenstate is the qq 
Fock component of the photon. The total 7*-proton 
cross section is easily calculated, if one knows the cross 
section <J qq (p) for scattering a qq-dipole of transverse 
size p off a proton, 
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The light-cone wavefunction ^> qq {a, p) describes the 
splitting of the virtual photon into the qq-pair and is 
calculable in perturbation theory, see e.g. [15]. Here, 
a is the longitudinal momentum fraction carried by 
the quark in Fig. 1. The dipole cross section is gov- 
erned by nonperturbative effects and cannot be cal- 
culated from first principles. We use the phenomeno- 
logical parameterization from [16] for our calculation 
of shadowing in DIS, which is fitted to HERA data 
on the proton structure function. Note that higher 
Fock-states of the photon, containing gluons, lead to 
an energy dependence of o qqi which we do not write 
out explicitly. 

The diffractive cross section can also be expressed 
in terms of a qq . Since the cross section for diffrac- 
tion is proportional to the square of the T-matrix ele- 
ment, |(7*|T|X}| 2 , the dipole cross section also enters 
squared, 
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where the brackets (...) indicate averaging over the 
light-cone wavefunctions like in Eq. (4). We point out 



that in order to reproduce the correct behavior of the 
diffractive cross section at large Mx, one has to in- 
clude at least the qqG Fock-state of the 7*. This cor- 
rection is, however, of minor importance in the region 
where shadowing data are available. 

If one attempts to calculate shadowing from Eq. (1) 
with help of Eq. (5), one faces the problem that the 
nuclear form factor, Eq. (2), depends on the mass Mx 
of the diffractively produced state, which is undefined 
in impact parameter representation. Only in the limit 
I c ^> Ra, where Ra is the nuclear radius, it is possible 
to resum the entire multiple scattering series in an 
cikonal- formula 

?qq{p), 
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The nuclear thickness function T(b) = dz pA(b, z) 
is the integral of nuclear density over longitudinal co- 
ordinate z and depends on the impact parameter b. 
The condition l c 3> Ra insures that the p does not 
vary during propagation through the nucleus (Lorentz 
time dilation) and one can apply the eikonal approxi- 
mation. 

The condition l c ^> Ra is however not fulfilled in 
experiment. For the case l c ~ Ra, one has to take the 
variation of p during propagation of the qq fluctuation 
through the nucleus into account, see Fig. 3. A widely 
used recipe is to replace M x — > Q 2 , so that l c — > 
1/ ' (2rriNXBj) and one can apply the double scattering 
approximation. This recipe was, however, disfavored 
during our investigation [17]. Moreover, there is no 
simple recipe to include a finite l c into higher order 
scattering terms. 




FIG. 3: Propagation of a gg-pair through a nucleus. 
Shown is the case of a finite coherence length, where 
the transverse motion is described by the Green function 
W (p 2 ,z 2 \f> 1 ,z 1 ). 

In [18] a Green function technique was developed 
that provides the correct quantum-mechanical treat- 
ment of a finite coherence length in all multiple scat- 
tering terms. Like in Eq. (1) the total cross section is 
represented in the form 
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is the minimal longitudinal momentum transfer when 
the photon splits into the qq dipole (m/ is the quark 
mass) . 

The shadowing term in Eq. (7) is illustrated in 
Fig. 3. At the point z\ the photon diffractively pro- 
duces the qq pair (7*^ — ► qqN) with a transverse 
separation p\. The pair propagates through the nu- 
cleus along arbitrarily curved trajectories, which are 
summed over, and arrives at the point Z2 with a sep- 
aration p 2 . The initial and the final separations are 
controlled by the light-cone wavefunction ^ q q(p, a). 
While passing the nucleus the qq pair interacts with 
bound nucleons via the cross section <J q q(p) which de- 
pends on the local separation p. The Green function 
W(p 2 , Z2\pi, Z\) describes the propagation of the pair 
from z\ to z 2 , see Eq. (9), including all multiple rcscat- 
terings and a finite coherence length. Note the diffrac- 
tion dissociation (DD) amplitude, 

Jdd{i* qq) = i*gg(pi,a)ov,g(/9i) (11) 

, in Eq. (9). At the position z 2 , the result of the prop- 
agation is again projected onto the diffraction dissoci- 
ation amplitude, Eq. (8). The Green function includes 
that part of the phase shift between the initial and the 
final photons which is due to transverse motion of the 
quarks, while the longitudinal motion is included in 
Eq. (9) via the exponential. 

The Green function W(p 2 , z 2 ; pi,z\) in Eq. (9) sat- 
isfies the two dimensional Schrodinger equation, 
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with the boundary condition W(p2, Z\ \p±, Z\) = 
^ 2 \p2 — Pi)- The Laplacian A(p 2 ) acts on the co- 
ordinate p 2 . The kinetic term A/[2i/a(l — a)] in this 
Schrodingcr equation takes care of the varying effec- 
tive mass of the qq pair and provides the proper phase 
shift. The role of time is played by the longitudinal 
coordinate z^. The imaginary part of the optical po- 
tential describes the rescattering. 

The Green function method contains the eikonal ap- 
proximation Eq. (6) and the Karmanov-Kondratyuk 
formula Eq. (1) as limiting cases. In order to ob- 
tain the eikonal approximation, one has to take the 
limit v — > oo. In this case, the kinetic energy term 
in Eq. (12) can be neglected and with q™ m — > one 
arrives after a short calculation at Eq. (6). 

One can also recover the Karmanov-Kondratyuk 
formula, when one neglects the imaginary potential 
in Eq. (12). Then W becomes the Green function of 
a free motion. 
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FIG. 4: The XBj dependence of shadowing in DIS for the 
structure function of tin relative to carbon. The data are 
from NMC [1]. The full curves are calculated in the Green- 
function approach, including the nonperturbative interac- 
tion between the q and the q. The dashed curve does not 
include this interaction. The figure is from [16]. 

Calculations in the Green function approach are 
compared to NMC data in Fig. 4. Note that the lead- 
ing twist contribution to shadowing is due to large 
dipole sizes, where nonperturbative effects, like an 
interaction between the q and the q, might become 
important. Therefore, the solid curve is calculated 
following [19], including such an interaction. This 
interaction modifies the light-cone wavefunction [19]. 
The dashed curve is calculated with the conventional, 
perturbative light-cone wavefunctions, but including a 
constituent quark mass. Both curves are in reasonable 



agreement with the data. 

Note that for these data, the coherence length is of 
order of the nuclear radius or smaller. Indeed, shad- 
owing vanishes around XBj ~ 0.1, because the coher- 
ence length becomes smaller than the mean internu- 
cleon spacing. Therefore, the eikonal approximation, 
Eq. (6), cannot be applied for the kinematics of NMC 
and a correct treatment of the coherence length be- 
comes crucial. We emphasize that the calculation in 
Fig. 4 does not contain any free parameters. Follow- 
ing the spirit of Glauber theory, all free parameters 
are adjusted to DIS off protons. 



3. NUCLEAR EFFECTS IN THE DRELL-YAN 
PROCESS 

In the second part of this paper, we discuss the ap- 
plication of the dipole approach to calculate nuclear 
effects in DY. The main features of the approach have 
already been introduced in the preceding section and 
the case of DY is almost completely analogous to DIS. 

The cross section for radiation of a virtual photon 
from a quark after scattering on a proton, can be writ- 
ten in factorized light-cone form [6, 7, 8], 



<Pp\^. q (a,p)\ 2 <T qq (ap), (13) 



da(qp — ► 7*Y) 
din a 

with the same a qq as in DIS. Again, the light-cone 
wavefunction ^> 1 t q {a,p) is calculable in perturbation 
theory [11, 20]. For DY, however, p is the photon- 
quark transverse separation and a is the fraction of 
the light-cone momentum of the initial quark taken 
away by the photon. We use the standard notation for 
the kinematical variables, x\ — X2 = xp, r = M 2 / s = 
X1X2, where xf is the Feynman variable, s is the center 
of mass energy squared of the colliding protons and M 
is the dilepton mass. 

The physical interpretation of (13) is similar to the 
DIS case. Long before the target, the projectile quark 
develops fluctuations which contain virtual photons. 
On interaction with the target, these fluctuations can 
be put on mass shell and the 7* is freed. In or- 
der to make this happen, the interaction must distin- 
guish between Fock states with and without 7*, i.e. 
these states have to interact differently. Since only the 
quark interacts in both Fock components, the differ- 
ence in the scattering amplitudes arises from the rela- 
tive displacement of the quark in the transverse plane 
when it radiates a photon. The j*q fluctuation has a 
center of gravity in the transverse plane which coin- 
cides with the impact parameter of the parent quark. 
One finds that the transverse separation between the 
quark and the center of gravity is ap, i.e. the argu- 
ment of a qq . The dipole cross section arises from the 
interference of the two graphs in Fig. 2 and should 



be interpreted as freeing cross section for the virtual 
photon. More discussion can be found in [20]. 

The transverse momentum distribution of DY pairs 
can also be expressed in terms of the dipole cross sec- 
tion [8]. The differential cross section is given by the 
Fourier integral 

d<r(qp — > 7*X) 
din ad 2 qr 

= J d2 Pi d2 P2 cxp[ig T • (pi - p 2 )] 

x *;» g (a,pi)* rg (a,p 2 ) (14) 

X 2 i a QQ( a Pl) + (J qq( a P2) - Vqq(a{pi - p 2 ))} . 

After integrating this expression over the transverse 
momentum qr of the photon, one obviously recovers 
Eq. (13). 

Nuclear effects arise in the DY process, because the 
different nucleons in the nucleus compete to free the 
7*. Necessary condition for this phenomenon is of 
course that the coherence length for DY, 
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exceeds the mean internucleon spacing of ~ 2 fm. At 
infinite coherence length, one can simply cikonalizc 
the dipole cross section, similar to the DIS case, cmp. 
Eq. (6). For finite coherence length, the Green func- 
tion technique can be developed for DY in complete 
analogy to DIS. We do not display the correspond- 
ing formula here, which can be found in [8, 16, 21]. 
We point out that even though DIS and DY appear 
very similar in the dipole approach, there are no (or 
only very few) diffractively produced DY pairs, be- 
cause Eq. (5) is not applicable to an exclusive channel 
(see discussion in [19]). 

Calculations of shadowing for DY in the Green func- 
tion approach are compared to E772 data [2] in Fig. 5. 
Like in the case of the NMC data in the previous sec- 
tion, the coherence length l c at E772 energy becomes 
smaller than the nuclear radius. Shadowing vanishes 
as X2 approaches 0.1, because the coherence length 
becomes smaller than the mean internucleon separa- 
tion. Again, it is therefore important to have a correct 
description of a finite l c . The eikonal approximation 
does not reproduce the vanishing shadowing toward 
x 2 — > 0.1, because it assumes an infinite coherence 
length. 

The curves in Fig. 5 arc somewhat different from 
the ones in [16, 21], because we used a different pa- 
rameterization of the dipole cross section [22]. Note 
that for heavy nuclei, energy loss [23] leads to an addi- 
tional suppression of the DY cross section. It is there- 
fore important to calculate shadowing in DY without 
introducing additional parameter for nuclear targets. 
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FIG. 5: Comparison between calculations in the Green 
function technique (solid curve) and E772 data at center 
of mass energy y/s = 38.8 GeV for shadowing in DY. The 
dashed curve shows the eikonal approximation, which is 
not valid at this energy, any more. 



Otherwise, shadowing and energy loss cannot be dis- 
entangled. 

Nuclear effects on the ^-differential cross section 
calculated at RHIC and LHC energy are shown in 
Fig. 6. At these high cm energies, the coherence length 
significantly exceeds the nuclear radius, justifying the 
eikonal approximation. However, in the case of such 
long l c , the quark has sufficient time to radiate ad- 
ditional gluons. Thus, higher Fock states have to be 
taken into account. This is accomplished by multiply- 
ing Oqq with the gluon shadowing ratio Rg , which was 
calculated in [19]. See [24] for a detailed discussion. 

The differential cross section is suppressed at small 




q x / GeV 

FIG. 6: Nuclear effects on the DY transverse momentum 
distribution at RHIC and LHC for dilepton mass M = 4.5 
GeV and Feynman xf = 0.5. 



transverse momentum qj- of the dilcpton, where large 
values of p dominate. This suppression vanishes at 
intermediate q? ~ 2 GeV. The Cronin enhancement 
that one could expect in this intermediate qr region 
[21] is suppressed due to gluon shadowing [24] . At very 
large transverse momentum nuclear effects vanish. 

4. SUMMARY 

In the target rest frame and at low XBj (#2), the 
cross sections for DIS and DY can both be expressed in 
terms of the same cross section for scattering a color- 
neutral gg-pair off a proton. The advantage of this 
formulation and the main motivation to develop this 
approach is the easy generalization to nuclear targets 
and the insight into the dynamical origin of nuclear 
effects, which appear as coherence effects due to mul- 
tiple scattering. 

The main nonperturbative input to all formulae is 
the dipole cross section, which cannot be calculated 
from first principles. Instead one uses a phenomeno- 
logical parameterization for this quantity, which is 
determined to fit DIS data from HERA. Nuclear ef- 
fects are then calculated without introducing any new 
parameters. A parameter free calculation of nuclear 
shadowing for DY is indispensable, if one aims at ex- 
tracting the energy loss of a fast quark propagating 
through nuclear from DY data [23] . Furthermore, the 
dipole formulation is the only known way to calcu- 
late the Cronin effect in a parameter free way. This 
approach can also be applied to calculate the Cronin 
effect in hadronic collisions [25]. Furthermore, a good 
description of NMC data on shadowing in DIS is 
achieved. 

Finally, we mention that the dipole approach can 
also be applied to a variety of other processes, e.g. 
heavy quark, quarkonium or vector meson production 
[26]. 
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